Chaos can act as a decoherence suppressor 
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We propose a strategy to suppress decoherence of a solid-state qubit coupled to non-Markovian noises by 
attaching the qubit to a chaotic setup with the broad power distribution in particular in the high-frequency do- 
main. Different from the existing decoherence control methods such as the usual dynamics decoupling control, 
high-frequency components of our control are generated by the chaotic setup driven by a low-frequency field, 
and the generation of complex optimized control pulses is not necessary. We apply the scheme to supercon- 
ducting quantum circuits and find that various noises in a wide frequency domain, including low-frequency 
1//, high-frequency Ohmic, sub-Ohmic, and super-Ohmic noises, can be efficiently suppressed by coupling the 
qubits to a Duffing oscillator as the chaotic setup. Significantly, the decoherence time of the qubit is prolonged 
approximately 100 times in magnitude. 
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Introduction. — Solid state quantum information process- 
ing 111] develops very rapidly in recent years. One of the basic 
features that makes quantum information unique is the quan- 
tum parallelism resulted from quantum coherence and en- 
tanglement. However, the inevitable interaction between the 
qubit and its environment leads to qubit-environment entan- 
glement that deteriorates quantum coherence of the qubit. In 
solid state systems, the decoherence process is mainly caused 
by the non-Markovian noises induced, e.g., by the two-level 
fluctuators in the substrate and the charge and flux noises in 
the circuits 0-0]. 

There have been numbers of proposals for suppressing 
non-Markovian noises in solid-state systems. Most of them 
suppress noises in a narrow frequency domain, e.g., low- 
frequency noises J2, [H. Among the proposed decoher- 
ence suppression strategies, the dynamical decoupling con- 
trol (DDC) [5] is relatively successful in suppressing non- 
Markovian noises in a broad frequency domain and has re- 
cently been demonstrated in the solid-state system experimen- 
tally [g]. The main idea of the DDC is to utilize high fre- 
quency pulses to flip states of the qubit rapidly, averaging out 
the qubit-environment coupling. The higher the frequency of 
the control pulse is, the better the decoherence suppression 
effects are. Efforts have been made to optimize the control 
pulses [7] in the DDC, however, the requirements of generat- 
ing extremely high frequency control pulses or complex op- 
timized pulses make it difficult to be realized in solid-state 
quantum information system, in particular in superconducting 
circuits. 

In this letter, we propose a method to extend the decoher- 



ence time of the qubit by coupling it to a chaotic setup fU • Al- 
though it is widely believed that the chaotic dynamics induces 
inherent decoherence 10-13], e.g., the quantum Loschmidt 
echo lH, we find surprisingly that the frequency shift of 
the qubit induced by the chaotic setup, which has not drawn 
enough attention in the literature, can help to suppress deco- 
herence of the qubit. 

The main merits of this method are: (1) the high frequency 
components, which contribute to the suppression of the non- 
Markovian noises, can be generated by the chaotic setup even 
driven by a low-frequency field; and (2) generating complex 
optimal control pulses is not necessary. 



Decoherence suppression by chaotic signals. - 
the following system-environment model ||4j] 
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where uj q (u>i), S z (a z ), and S± (a 1 ^) are the angular fre- 
quency, the z-axis Pauli operator, and the ladder operators of 
the qubit (the i-th two-level fluctuator in the environment); gi 
is the coupling constant between the qubit and the i-th two- 
level fluctuator; and S q (t) = S q is an angular frequency shift 
induced by chaotic signals. If the initial state of the system is 
separable, p(0) = pso <8> Pbq, we can reduce the dynamical 
equation of the total system by tracing out the degrees of free- 
dom of the bath. The influence of the chaotic signal S q (t) falls 
into two aspects: (1) S q (t) affects the angular frequency shift 
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of the qubit induced by the bath 



Hamiltonian 



Aw, = 



Im 
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and (2), more importantly, it modifies the bath-induced deco- 
herence rate of the qubit which can be written as: 
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where J (u) = £\ g 2 <5 (w 



is the spectral density of the 



bath. Here, since the frequencies of the fiuctuators distribute 
in a finite domain, Aw g and T q are restricted to be integrated 
in the finite frequency domain [uj c i, ^>c2\- We demonstrate our 
results using the zero-temperature bath. 

We now come to show how the damping rate can be reduced 
by the frequency shift 8 q {t), using the function 8 q (t) as a lin- 
ear combination of sinusoidal signals with small amplitudes 
and high frequencies, i.e., 6 q (t) = J2 a Ada cos (u!d a t + <p a ). 
Here ujda should satisfy the conditions: ujda ^ |^c 2 — w g|> 
\co q — uj Ci |, and Ada da *C 1. Using the Fourier-Bessel se- 
ries identity H]: e" siny = J n 0) e iny with J n (x) as 



the n-th Bessel function of the first kind and the approxima- 
tion Jq {%) 1=3 1 — x 2 /4 for x -C 1, we have 
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where the correction factor F = cxp ^— 7r JJ°^ 

£<5 9 M = S Q -^da* - w ) / 2?r is the P ower spectrum 
density of the signal S q (t); uj- — (ui q — uf) /2; and u c d is 
the lower bound of the frequency of S q (£) such that ui c d 
\cu C2 — u> q \, \oj q — u> Cl \. Here, we omit the higher-order 
Bessel function terms because 1/ (i (uj q — ui) + noJda) 
1/ (i (uj q - u)) and J n {A da /udu) < Jo (A da /u da ) under 
the conditions. The analysis shows that 6 q (t) induces a cor- 
rection factor F for the environment-induced frequency-shift 
Aui q and damping rate T q , i.e., 



Auj q = F 



r„(t) = f 



^JMti-cosK-^f] = 

2 (Wq - LO) 

d JlJ (w)sin (cjg--^)j = ^ 



grO i 



W 



(/< ) , 



where Aw 9 o and T q Q are the frequency-shift and damping 
rate when 8 q (t) = 0. The correction factor F may become 
extremely small when 6 q (t) is a chaotic signal which has 
a broadband frequency spectrum in particular in the high- 
frequency domain, such that the decay rate and frequency shift 
can be suppressed by a chaotic signal. 

Generation of chaotic signals and suppressing 1/f 
noises. — To show the validity of our method, as an exam- 
ple, we show how to suppress the 1/f noises of a qubit with 
free Hamiltonian H q = u) q S z by coupling it to a driven Duff- 
ing oscillator which is used to generate chaotic signals, with 



Duf = L0 o a<a- - [a 



(2) 



and damping rate 7, where a and fit are the annihilation and 
creation operators of the nonlinear Duffing oscillator; w /2ir 
is the frequency of the fundamental mode; A is the nonlinear 
constant; and / (t) = Jo cos (uidt) denotes the classical driv- 
ing field with amplitude Iq and frequency Ud/ln. We employ 
the interaction between the qubit and the Duffing oscillator, 
Hi = g qo S z tJa (g qo - coupling strength), which can be ob- 
tained, e.g., by the Jaynes-Cummings model under the large 
detuning regime 11 1 511 . 

By tracing out the degrees of freedom of the oscillator ini- 
tially in a coherent state |a), we find that the interaction be- 
tween the qubit and the Duffing oscillator introduces an ad- 
ditional factor for the non-diagonal entries of the state of the 
qubit @]: 



foi (t) 
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There are two aspects of the factor /01 if) = e 



E a (t)+»e„(t) 



The phase shift Q q (t) w J Q S q (t')dt', induced by foi (t), is 
related to S q (t) in Eq. (Q]), which can be used to suppress 
the decoherence of the qubit. However, the amplitude square 
M{t) = l/m (t)\ 2 = e 2S « (t) , i.e., the quantum Loschmidt 
echo iTToL Flltl . leads to additional decoherence effects of the 
qubit. Such decoherence effects have been well studied in the 
literature for regular and chaotic dynamics 10-13]. Our de- 
coherence suppression strategy is valid when the decoherence 
suppression induced by S q (t) is predominant in comparison 
with the opposite decoherence acceleration process. 

We now come to show numerical results, using system pa- 
rameters: 



K, g qo , 1, A) = (u q , 0.03w g , 0.05w 9 , 0.25o; g ) . 



(4) 



The bath has a 1// noise spectrum (see, e.g., Ref. |2j,|3[]) with 
J (oS) = A/uj and A/oj q = 0.1. The evolution of the co- 
herence C xy = (S x ) 2 + (S y ) 2 of the qubit and the spectrum 
analysis of the angular frequency shift 5 q (t) are presented in 
Fig-Ill As shown in Fig.QJb), (c), if we tune the amplitude I 
of the sinusoidal driving field I (t) such that Io/uJq = 5 and 
30, the signals S q (t) exhibit periodic and chaotic behaviors. 
As shown in Fig. 0I a X m the periodic regime, the decoher- 
ence of the qubit is almost unaffected by the Duffing oscilla- 
tor. The trajectory in the periodic case (green curve with plus 
signs) coincides with that of natural decoherence (black trian- 
gle curve), as in Fig.[TJa). In the chaotic regime, the decoher- 
ence of the qubit is efficiently slowed down (see the blue solid 
curve in Fig. [TJa) representing the trajectory in the chaotic 
case). This demonstrates that, with the increase of the distri- 
bution of the spectral energy in the high-frequency domain, 
the decoherence effects are suppressed as explained in the last 
section. 
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FIG. 1: (color online) Decoherence suppression by the auxiliary 
chaotic setup, (a) the evolution of the coherence C xy = (&) 2 + 
(Sy) 2 of the state of the qubit, where the red asterisk curve and the 
black triangle curve represent the ideal trajectory without any de- 
coherence and the trajectory under natural decoherence and without 
corrections; and the green curve with plus signs and the blue solid 
curve denote the trajectories with Io/u) q = 5 and 30. With these 
parameters, the dynamics of the Duffing oscillator exhibits periodic 
and chaotic behaviors, r = 2n/uj q is a normalized time scale, (b) 
and (c) are the energy spectra of 5 q (t) with Io/uiq ~ 5 (the peri- 
odic case) and 30 (the chaotic case). The energy spectrum Ss q (w) 
is in unit of decibel (dB). (d) the normalized decoherence rates V /uj q 
versus the normalized driving strength Io/uj q . 

We compare in Fig. |TJd) the average natural decoherence 
rate r 9 o and modified decoherence rate T qm of the qubit ver- 
sus different strengths I of the driving field. Figure Q~fd) 
shows that, when the strength I Q of the driving field increases, 
the decoherence process is efficiently slowed down. It is 
interesting to note that there seems to exist a phase transi- 
tion around Io/u>q = 20, i.e., a sudden change of the mod- 
ified decoherence rate T qm (see the black solid curve with 
plus signs). It is noticable that, around this point, the dy- 
namics of the Duffing oscillator enters the chaotic regime 
which is indicated by a positive Lyapunov exponent (see the 
green dash-dotted curve in Fig.QJd)). The modified decoher- 
ence rate T qm changes dramatically in the parameter regime 
Io/ijJq € [20, 35] which is the soft-chaos regime of the Duff- 
ing oscillator. When the dynamics of the Duffing oscillator 
enters the hard-chaos regime at Io/u> q rj 35, the modified de- 
coherence rate T qm is stabilized at a value much smaller than 
the natural decoherence rate r g o- The simulation results show 
that the decoherence of the qubit is efficiently suppressed by 
our proposal, even if there exists an additional decoherence in- 
troduced by the auxiliary chaotic setup. Further calculations 
show that the modified decoherence rate T qm in the chaotic 
regime is roughly 100 times smaller than the unmodified de- 
coherence rate T q o, meaning that the decoherence time of the 
qubit can be prolonged 100 times. 

Experimental feasibility in superconducting circuits. — Our 
general study can be demonstrated using the solid state 



quantum devices, e.g., the superconducting qubit system, as 
sketched in Fig. [2] It is similar to the widely used qubit 
readout circuit II 1 6fl . but works in a quite different parame- 
ter regime. In this superconducting circuit, a single Cooper 
pair box (SCB) is coupled to a dc-SQUID consisting of two 
Josephson junctions with capacitances Cj and Josephson en- 
ergies Ej and a paralleled current source. The SCB is com- 
posed of two Josephson junctions with capacitances Gj and 
Josephson energies Ej, The difference between the circuit in 
Fig. |2]and the readout circuit in Ref. II 1611 is that the rf -biased 
Josephson junction is replaced by a dc-SQUID - the chaotic 
setup. 




FIG. 2: (color online) Schematic diagram of the decoherence sup- 
pression superconducting circuit in which a SCB is coupled to a 
current-biased dc-SQUID. 

The Hamiltonian of the circuit shown in Fig.|2]can be writ- 
ten as: 

H = E c {n~ n g f - 2Ej cos - cos 9 + E c n 2 

— IE j cos cos <f> — <poI e <p, (5) 

where Ec = 2e 2 / [C g + 2Cj) is the charging energy of 
SCB with C g as the gate capacitance; n g = —C g V g /2e is 
the reduced charge number, in unit of the Cooper pairs, with 
V g as the gate voltage; n is the number of Cooper pairs on 
the island electrode of SCB with 9 as its conjugate opera- 
tor; Eq = e 2 jCj is the charging energy of the dc-SQUID; 
n is the charge operator of the dc-SQUID and <p the conjugate 
operator; and <f) e and I e are the external flux threading the 
loop of the dc-SQUID and the external bias current of the dc- 
SQUID. Here, we consider a zero external flux threading the 
loop of the coupled SCB-dc SQUID system. In this case, the 
phase drop across the SCB is equal to the phase drop across 
the dc-SQUID <f>. We further introduce the ac gate voltage 
Vg = V g o cos (uigt) with amplitude V g o and angular frequency 
id g . With the condition that C g V g oEc /2e <C u> q = Ej — ui g , 
the SCB works near the optimal point, and we only need 
to worry about the relaxation of the SCB. By expanding the 
Hamiltonian of the SCB in the Hilbert space of its two low- 
est states and leaving the lowest nonlinear terms of (j> in 
the rotating frame, we can obtain the effective Hamiltonian 
H c s = H q + Huui + Hi discussed in the foregoing sec- 
tion. This dc-SQUID, acting as the auxiliary Duffing oscil- 
lator, can be used to suppress low frequency 1// noises of 
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Type of noises 


Frequency domain 


r g o 




Ti =T 2 


1/f noise 
Ohmic 
Sub-Ohmic 
Super-Ohmic 


[10 kHz, 1 MHz] 

[2^/3,3^/2] 
[2cV3,3(V2] 
[2^/3,3^/2] 


0.58 MHz 
0.35 MHz 
0.35 MHz 
0.36 MHz 


5.4 kHz 
5.4 kHz 
5.4 kHz 
5.4 kHz 


187 us 
187 fis 
187 fis 
187 /is 



TABLE I: Decoherence suppression against various noises for exper- 
imentally accessible parameters: Ej/2tt = 5 GHz, u g /2n = 4.999 
GHz, E c /2tt = 0.188 MHz, and Ej cos if /2n = 12.032 MHz. 



60904034, 60836001, 60635040. T. J. Tarn would also like to 
acknowledge partial support from the U. S. Army Research 
Office under Grant W91 1NF-04- 1-0386. L. A. Wu has been 
supported by the Ikerbasque Foundation Start-up, the Basque 
Government (grant IT472-10) and the Spanish MEC (Project 
No. FIS2009-12773-C02-02). 



the qubit. Using the experimentally accessible parameters as 
shown in the caption of table|I] we show the decoherence sup- 
pression effects for low frequency 1/f, high frequency Ohmic 
(J(cj) = we-"/ 5 ""), sub-Ohmic (J (w) = uP-^e'"^), 
and super-Ohmic (J (uj) = ui 2 e~ ul / buJq ) noises. All simula- 
tions are summarized in table [I] It is found that our method 
works equally well for different types of noises. The numeri- 
cal simulations manifest that our strategy is independent of the 
sources and frequency domains of the noises. The final modi- 
fied decoherence rates for these different noises are almost the 
same because the decoherence effects induced by the environ- 
mental noises are all greatly suppressed, and thus the modi- 
fied decoherence rates of the qubit are mainly caused by the 
auxiliary chaotic setup, i.e., the dc-SQUID. It is also shown 
in table U that the modified decoherence rate T qm /2Ti of the 
qubit can be reduced to 5 kHz. This low decoherence rate cor- 
responds to a long decoherence time T\ = T2 ~ 200 /is. The 
magnitude is one-order longer than the decoherence time of 
the superconducting qubits realized in experiments (see, e.g., 
Ref. H). 

Conclusion. — In conclusion, we propose a strategy to pro- 
long the decoherence time of a qubit by coupling it to a chaotic 
setup. The broad power distribution of the auxiliary chaotic 
setup in particular in the high-frequency domain helps us to 
suppress various non-Markovian noises, e.g., low-frequency 
1/f noise, high-frequency Ohmic, sub-Ohmic, and super- 
Ohmic noises, and thus freeze the state of the qubit even if we 
consider the additional decoherence induced by the chaotic 
setup. We find that the decoherence time of the qubit can be 
efficiently prolonged approximately 100 times in magnitude. 
We believe that our strategy is feasible, in particular for a cou- 
pled SCB-SQUID system, and also gives a new perspective 
for the reversibility and irreversibility induced by nonlinear- 
ity. 

J. Zhang would like to thank Dr. H. T. Tan and Dr. M. H. 
Wu for helpful discussions. We acknowledge financial 
support from the National Natural Science Foundation of 
China under Grant Nos. 60704017, 10975080, 61025022, 



* Electronic address: jing-zhang@mail.tsinghua.edu.cn 
[1] J. Q. You and F. Nori, Physics Today 58 (11), 42 (2005); 

Y. Makhlin et al, Rev. Mod. Phys. 73, 357 (2001); J. Clarke 

and F. K. Wilhelm, Nature 453, 1031 (2008). 
[2] F. Yoshihara et al, Phys. Rev. Lett. 97, 167001 (2006); R. C. 

Bialczak et al, ibid. 99, 187006 (2007); K. Kakuyanagi et al, 

ibid. 98, 047004 (2007); L. Faoro and L. Viola, ibid. 92, 1 17905 

(2004); E. Collin et al, ibid. 93, 157005 (2004); L. Tian, ibid. 

98, 153602 (2007); R. H. Koch et al, ibid. 98, 267003 (2007). 
[3] G. Ithier et al, Phys. Rev. B 72, 134519 (2005); J. M. Martinis 

et al, ibid. 67, 094510 (2003); K. Shiokawa and D. A. Lidar, 

Phys. Rev. A 69, 030302(R) (2004); G. Falci et al, ibid. 70, 

040101 (R) (2004). 
[4] J. Bergli et al, New J. Phys. 11, 025002 (2009). 
[5] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998); L. Viola 

et al, Phys. Rev. Lett. 82, 2417 (1999). 
[6] J. F. Du et al, Nature 461, 1265 (2009). 

[7] G. S. Agarwal, Phys. Rev. A 61, 013809 (1999); L.-A. Wu and 
D. A. Lidar, Phys. Rev. Lett. 88, 207902 (2002); A. G. Kofman 
and G. Kurizki, Phys. Rev. Lett. 93, 130406 (2004); G. S. Unrig, 
Phys. Rev. Lett. 98, 100504 (2007); W. Yang and R.-B. Liu, 
Phys. Rev. Lett. 101, 180403 (2008); W. Yao et al, Phys. Rev. 
Lett. 98, 077602 (2007). 
[8] E. N. Pozzo etal, Phys. Rev. B 77, 024518 (2008); E. N. Pozzo 
and D. Dominguez, Phys. Rev. Lett. 98, 057006 (2007); T. Mao 
and Y. Yu , Phys. Rev. E 81, 016212 (2010); N. Lambert et al, 
Phys. Rev. B 80, 165308 (2009); M. J. Everitt etal, Phys. Lett. 
A 374, 2809 (2010). 
[9] W.-G. Wang et al, Phys. Rev. A 77, 012108 (2008). 

[10] A. Peres, Phys. Rev. A 30, 1610 (1984); A. Peres, Quantum 
Theory: Concepts and Methods (Kluwer, Dordrecht, 1995). 

[11] R. A. Jalabert et al, Phys. Rev. Lett. 86, 2490 (2001); F. M. 
Cucchietti etal, Phys. Rev. Lett. 91, 210403 (2003); G. Veble 
et al, Phys. Rev. Lett. 92, 034101 (2004); H. T. Quan et al, 
Phys. Rev. Lett. 96, 140604 (2006). 

[12] J. Vanicek and E. J. Heller, Phys. Rev. E 68, 056208 (2003). 

[13] Ph. Jacquod et al, Phys. Rev. E 64, 055203(R) (2001). 

[14] L. Zhou et al, Phys. Rev. A 80, 062109 (2009). 

[15] A. Blais et al, Phys. Rev. A 69, 062320 (2004). 

[16] D. Vion et al, Science 296, 886 (2002); I. Siddiqi et al, Phys. 
Rev. Lett. 93, 207002 (2004); Phys. Rev. B 73, 054510 (2006). 



